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The spin-momentum locked nature of the robust surface states of three dimensional topological insulators (3D 
TI) make them promising candidates for spintronics applications. Surface potentials which respect time reversal 
symmetry can exist at the surface between a 3D TI and the trivial vacuum. These potentials can distort the 
spin texture of the surface states while retaining their gapless nature. In this work, the effect of all such surface 
potentials on the spin textures is studied. Since, a tunnel magnetoresistance signal carries the information of the 
spin texture, it is proposed that spin-polarized tunneling of electrons to a 3D TI surface can be used to uniquely 
identify the surface potentials and quantitatively characterize them. 


I. INTRODUCTION 

Topological insulators are exotic materials with an insu¬ 
lating bulk characterized by a gapped energy spectrum and 
topologically protected conducting surface states character¬ 
ized by massless Dirac fermionsThe surface states of 
strong three dimensional topological insulators (3D TI)^'^ 
are robust against disorder (which respect time reversal) in¬ 
duced backscattering as they are topologically protected by 
time reversal symmetry®’^°’*'’'^“*^. These surface states are 
special with regard to them being perfectly spin-momentum 
locked. This feature of the surface states have been confirmed 
experimentally using spin resolved angle resolved photoemis¬ 
sion spectroscopy (ARPES)*^’'®'**. Study of disorder induced 
quasi-particle interference patterns using a spin polarized 
scanning tunneling microscope (SP-STM)^®“^^ have also been 
used to probe this nature of surface states. Experiments with 
multiple ferromagnetic contacts have also been suggested^^’^^ 
and they do show indications of spin-momentum locking^"^^^®. 
More recently, there have also been proposals^^’^^ to recon¬ 
struct their spin texture of the surface states using spin po¬ 
larized injection of currents from a source like an SP-STM 
and extracting an effective tunnel magnetoresistance^® within 
a multi-terminal electrical transport setup. 

The spin-momentum locked nature of the surface states 
provide a lot of scope for manipulating and directing spin 
polarized current in a controlled fashion on the surfaces of 
these materials. This along with robustness of these sur¬ 
face states make them the ideal materials for applications to 
spintronics^®’^'. It then becomes extremely important to un¬ 
derstand any other effect that can potentially distort the spin 
texture. Indeed, there can exist potentials localized on the sur¬ 
faces of the 3D TIs which could affect the spin texture. It 
is then also desirable to have a protocol to characterize these 
surface potentials on a given sample. 

In principle, both time reversal symmetry breaking po¬ 
tential and potentials which respects time reversal could be 
present on the surface of a topological insulators. Observa¬ 
tion of robust surface state in experiments strongly indicate 
the absence of former but the presence of the later can not be 
ruled out in general. Hence studying the effect of such po¬ 
tentials can be of importance. In a recent work by Zhang et 
al. (Ref.32) an exhaustive list of such surface potentials was 
presented. Zhang et al. studied the effect of these potentials 


to explain doping related effect observed in experiments on 
3D TI surface states. In this article we will study the effect of 
such potentials on the spin texture of the surface states. 

Electrons on a planar surface of a 3D TI can be described in 
terms of two independent SU(2) degrees of freedom. For the 
surfaces perpendicular to the crystal growth are the electron 
spin and the orbital pseudospin. It turns out that these sur¬ 
face potentials can be classified into two categories, one that 
couples only to the orbital pseudospin of the electrons and the 
other which couples only to the spin. The second category of 
surface potentials which affect the spin textures is the focus of 
this work. 

It turns out that these surface potentials have non-trivial ef¬ 
fects on the surface states which lead to novel spin textures. 
For instance, some of the surface potentials tilt the spin out 
of the plane of the momentum leading to an imperfect spin- 
momentum locking. At some particular values of the sur¬ 
face potentials, the spin texture is distorted in such a way, 
that it looks as if it belongs to a Dresselhaus spin-orbit cou¬ 
pling even though the spin-orbit coupling term in the Hamil¬ 
tonian governing the surface states is of the Rashba-type®^. 
Another interesting effect of the surface potentials found was, 
the spin texture gets twisted in the momentum-plane such that 
the spin-momentum locking stays perfect, albeit with a new 
spin-momentum locking angle. In this work, we study the ef¬ 
fect of the surface potentials on the spin textures and provide 
analytical expressions for the spin textures in the presence of 
these potentials. 

The central message of this work is demonstrating how 
novel spin-textures can be obtained on the surface states of 3D 
TIs if the surface potentials can be manipulated in a controlled 
fashion. Further, using the protocol developed in Ref.27 we 
describe how these new spin textures could be identified using 
simple quantum transport experiments. In fact, our proposal 
also provides a direct route to the quantitative identification of 
the surface potentials present on a sample. 

The rest of the paper is organized as follows. In Sec. II, 
the effect of the surface potentials on the spin textures of the 
surface states is described. Sec. Ill briefly sketches the con¬ 
nection between spin-polarized tunneling current and the spin 
texture. Sec.IV demonstrates how spin-polarized tunneling 
can be used to extract information about the surface potentials 
via the spin textures and finally Sec.V discusses experimental 
feasibilities and summarizes the results. 
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II. EFFECT OF SURFACE POTENTIALS ON SPIN 
TEXTURE OF SURFACE STATES 

Electrons in 3D TI materials like Bi 2 Se 3 , Bi 2 Te 3 etc. have 
two SU(2) degrees of freedom, the orbital (r) and the spin 
(cr). The low energy physics of these materials is captured by 
the Hamiltonian 

T-Lbuik = ^ 2 ®{-mQTz+VzkzTy)+v\\{kyax-kx(Ty)®Tx ( 1 ) 

where it is assumed that the crystal growth axis of the sam¬ 
ple is along the z-direction. The boundaries of these crystals 
host topologically protected surface states as the mass under¬ 
goes as inversion in sign across the boundary from a 3D Tl 
to the trivial vacuum. These surface states are evanescent on 
both the sides of the boundary as dictated by the Jackiw-Rebbi 
solution^^ for the one-dimensional quantum mechanical prob¬ 
lem along a direction perpendicular to the surface. To ob¬ 
tain these surface states, one effectively needs to calculate the 
evanescent states at the boundary between the TI and the vac¬ 
uum at the z = 0 plane, using the Hamiltonian in Eq.l using 
a positive mass toq = to on the TI side (say z < 0) and a 
negative mass TOq = —M on the vacuum side (z > 0). Since, 
the vacuum is a trivial insulator, the limit of M —cx) is taken. 

The surface states so obtained are described by a massless 
Dirac Hamiltonian of the form 

'kisurf — t^ll ky (7^kx\ (2) 

whose eigenstates are perfectly spin-momentum locked. The 
Jackiw-Rebbi solution also says that the surface state on the 
z = 0 surface is an eigenstate of the operator with eigen¬ 
value +1^^. The continuity of the solution at z = 0 ensures 
that the spinor at z = 0 “, also has the same structure as that 
of the TI surface except that the decay length scale of these 
states on the two sides are different as dictated by the Dirac 
mass. Hereafter we denote this spinor as \tpo) and the spin 
texture evaluated from it turns out to be 

(cr) = (sin0k,-cos6»k,O). (3) 

where 0 k = tan.~^{ky/kx)- 

However, the presence of time reversal symmetric surface 
potentials can significantly affect the nature of the these sur¬ 
face states and consequently the spin texture. Note that, since 
the surface potentials respect time reversal symmetry, they do 
not destroy the surface states. 

It is important to stress that the surface potential cannot 
simply be incorporated into the surface Hamiltonian in Eq.2. 
In fact, to study the effect of the surface potentials on the 
surface states, it is required to re-derive the surface states by 
solving the Jackiw-Rebbi problem across the mass inversion 
boundary, but crucially now in the presence of the surface po¬ 
tentials. Hence, the effective problem now, is to obtain the 
evanescent modes at z = 0 of a system described by the 
Hamiltonian 

ki — Hbulk + Hpot, (4) 

where 

T-Lpot = (T®Ty5{z), ( 5 ) 


with A = (2uz/to)A and A = (A^;, Ay, A^). 

The spinor on the trivial vacuum side stays the same as \iPq) 
because the limit of M ^ -foo is anyway taken. However, 
the Dirac Hamiltonian being linear in the momentum opera¬ 
tor ( first order spatial derivative), the surface potentials mod¬ 
eled using the Dirac-delta functions (Eq.5) lead to discontinu¬ 
ities in the spinor across the surface between the TI and the 
trivial vacuum. As a result, the surface state spinors are not 
described anymore by |V'o)^ rather, they are described by a 
different spinor, denoted by |V’a)- It can be shown that the 
spinor \4’a) can be obtained from \tpQ) via a matrix A4 (see 
Appendix for the derivation) as 

Ii/'a) = AJIV'o), ( 6 ) 

which appropriately implements the topological boundary 
condition^^ and hence gives the discontinuities. The form of 
A4 is given by 

M = |\.|2 - 2imA • (T (g) I 2 ] . (7) 

A closer inspection of Eq.7 reveals that it can be written as 
a direct product of a rotation matrix in the a sector, and an 
identity matrix in the r sector as 

M = a 1 ,; 0 ,, = cos-' ( ”) ) l^j) ) ( 8 ) 

Hence the effect of the surface potentials can be geometrically 
interpreted as a rotation of the spin of the electron by an angle 
20 a about the axis A, while not affecting the orbital pseu¬ 
dospin. The identification of AJ as a rotation matrix in the 
a sector allows for a broad classification of the surface poten¬ 
tials into two classes, one which rotates the spin textures about 
an axis which lies in the plane of the surface i.e. 7 ^ 0 , 

Az = 0 and the other which rotates the spin textures about an 
axis perpendicular to the surface i.e. Ax(y) = 0, Az 7 ^ 0. The 
surface potentials of the first class can actually tilt the spin tex¬ 
ture out of the plane of the surface. The maximum angle by 
which it can rotate the spin is tt (| A| = to) when the spin tex¬ 
ture is back to being a planar one again. The important feature 
to note here is that, after a rotation by tt the spin-texture looks 
as if it comes from a Dresselhaus spin-orbit coupled Hamil¬ 
tonian. This extrema is indeed very special as it changes the 
winding number of the spin as one moves along the Eermi 
surface from -|-1 to — 1. The other class of surface potential is 
fundamentally different as it always rotates the spin texture in 
the plane of the surface and hence can never change the wind¬ 
ing number of the spin texture along the Eermi surface. It is 
also interesting to note that in both the limits of | A|/to ^ 1 
and IAI /to ^ 1 the angle 0 a —> 0 which points to the fact 
that there is an intermediate window in the magnitude of the 
surface potentials where it has a significant effect on the spin 
textures. It is important to note that since the surface poten¬ 
tials imposed are in the form of a direct product between the 
a and r spaces, the resultant M matrix is also of the same 
structure. Consequently the direct product form of the spinors 
(in the r and the a sector) for the surface states are retained. 
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The spin texture in the presence of the surface potentials can (cr^) = cos 0k + Zg sin 0k, (9c) 

be obtained by taking an expectation value of the cr operator 
as (V'aIo'IV'a)- The distorted spin textures turn out to be of 
the form 


(T^) = A’c cos 0 k + A’s sin 0 k, 

(9a) 


= yc cos 0 k - 1 - sin 0 k, 

(9b) 

where the coefficients of cos 0k and sin 0k are functions of A 
given by 


-T, = - 


4(Ai(A^2 + - 1) + 2A^A; + A'3) ^ 


(1 


,n 2'|2 


8(a;2 + a;2 ) 

(1+|A'|2)2 


( 10 a) 


8 (AL^ + A'/) 

(1 + |A'|2)2 ’ 


3^. 


4(a;(a;' + a;" -1) + 2a^a; + a'/) 

(1 + |A'|2)2 


(10b) 


z. = 


4(a;3 - 2a;a; + a;(-i + a;2 + a;^)) 

(1 + |A'|2)2 


where the A's are dimensionless numbers given by A^/m. 

The effect of the different surface potentials can be under¬ 
stood by looking at Eq.lO for the different cases separately. 
Firstly, in the absence of any surface potential (A = 0), Eq.9 
and Eq.lO reproduce the pristine spin texture mentioned in 
Eq.3. However, the presence of the surface potentials indeed 
leads to distorted spin textures. 

By setting A^ ^ 0 and A^ = 0 = A^ in Eq.lO the effect 
of the surface potential of the form can be understood. 

It turns out that such a potential leaves the spin texture along 
the cc-direction ((cr^)) undisturbed as A), = 0 and A’s = 1 
in this case. On the contrary, such a potential reduces the 
magnitude of the spin texture along the y-direction and also 
tilts the spin out of the x-y plane which is evident from (cr^) 
being non-zero. This is exactly the effect that a rotation of the 
spin texture about the x axis produces. The effect of surface 
potentials of the form Aya^ and A^cr^ can also be understood 
similarly from the interpretation of the rotation of the spin 
texture as mentioned in Eq. 8 . 


III. SPIN POLARIZED TUNNELING 

In this section, the calculation of spin-polarized tunneling 
current on the surface of a 3D TI is outlined and shown that 
how it carries information of the spin texture on the Fermi 
surface of the 3D TI surface states. In the second quantized 
language, the ground states of a 3D TI and a spin-polarized 
electron source can be written as Y[k<kF and 
Y[k<kF V’S'^k respectively where the ipTi(S) are the spinors 
describing the internal degrees of freedom and are 

usual fermionic creation operators for the 3D TI surface (spin 
polarized source). For a particular momentum mode, the spin 


Z. = 


4(a;3 + 2a;a; + a;(-i + a'^ + a'^)) 

■(1 + | A '| 2)2 


(lOc) 


of the 3D TI surface state and the spin polarized source can 
respectively be calculated as S 7 ’/(k) = (' 0 T/(k)|cr|'!/; 7 '/(k)) 
and S 5 = {'ll} s\(t\' 4> s) ■ Note that the spin of the spin-polarized 
source has no momentum dependence. 

The tunneling Hamiltonian between such a source and the 
3D TI surface can be written as 

Utunn = J (T'^(r = 0 ) 4 'T 7 (r = 0) + h.c). (11) 

In the momentum space Titunn looks like 

'^tunn — J H' h.c), (12) 

where J is the tunneling amplitude, is the overlap of the 
spinors of the 3D TI surface state and the spin polarized 
source. The overlap of the spinors has a clear geometrical 
interpretation as 


|zk|" = l(l + ST/(k)-Ss) (13) 

where the ST/(k) and S 5 are unit vectors representing the 
spin expectation values of the 3D TI surface state and the spin 
polarized source. Note that this is the well known tunnel mag¬ 
netoresistance response form for transport between two spin- 
polarized materials 2 ®. 

The total current turns out to be a sum of the momentum- 
resolved current values over the available momentum modes 
in the bias window2^’2^. This fact is extremely crucial as the 
decomposition of the total current into momentum-resolved 
currents enables us to divide the current that flows in each 
half of the TI surface to each of the contacts shown in Fig. 1(c) 
by simply summing the current over the momentum modes 
which go towards each contact from the spatial point r = 0 . 
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Presence of time reversal symmetry ensures that the net mag¬ 
netization of the entire Fermi surface is identically zero, hence 
a current response summed over the entire Fermi surface can¬ 
not give any information about the spin texture. However, 
finite segments of the Fermi surface do have non-zero mag¬ 
netization leading to a spin dependent response to tunneling 
which carry information about the spin texture, hence being 
able to separate the current response from segments of the 
Fermi surface is necessary to extract any information of the 
spin texture. 

The expression for the momentum resolved current is given 

by27,28 

/(k)=®^(l + ST/(k).Ss)x 
a 

(nF(i?k,/2 t7, Tt/) — ^s,Ts)), (14) 

where is the density of the states of the spin-polarized 
source. The momentum resolve current allows us to define 
a dimensionless current symmetry (denoted by A/) which fil¬ 
ters out the information of the spin texture. For a given config¬ 
uration of the contacts determined by the angle 7 (see Fig.lc), 
A/ is defined as 

C 7+7r 7+27r\ oo 

/ - / UOk JdkkKk) 

" 7 *'2 -^- (15) 

27r 00 

/ d6»k J dkk /(k) 

0 0 


The way AI is defined makes it couple to the magnetization 
of half of the Fermi surface. Time reversal symmetry ensures 
that the magnetization of one half of the Fermi surface is al¬ 
ways equal and opposite to that of the other half. The isotropy 
of the Fermi surface of the surface states of 3D TI tell us that 
the spin texture depends just on 9^ and the Fermi energy just 
on |k|. As a result, the current asymmetry AI takes a simple 
form given by 

2 rJ+Tk 2 

AI = -Ss- / d9i^STi{9k) = —Ss-STi.haif{"f), (16) 

TT TT ' 

where STi.haifil) denotes the magnetization of half of the 
Fermi surface between the azimuthal angles 7 and 7 -f tt. 

To illustrate how the dimensionless current asymme¬ 
try defined in Eq.l6 reconstructs the spin texture^^, the 
example of the spin texture of a pristine surface state 
(Eq.3) is shown in Eig.l. The spin of the spin-polarized 
source can be parametrized using its polar angles Sg = 
(sin 0s cos (/)s, sin 0s sincos 0s). As mentioned in 
Ref.27, the spin-momentum locking angle, 0^ for a planar 
spin texture can be obtained by setting 0s = 7r/2 and ob¬ 
serving the angle ’jmax at which AI shows an extrema as 
Imax = — ^L+7r/2. Note that the same magnitude of 

AI at ‘Jmax for (j>s — 0 and (pg = 7r/2 in Eig.lc is a signa¬ 
ture of the perfect spin-momentum locking. Also, the signs of 
the extrema tell us about the chirality and type of spin-orbit 
coupling, Rashba or Dresselhaus. In order to detect and re¬ 
construct an out of plane spin polarization in the Tl surface 
state, one simply needs to set 0s = 0. 



FIG. 1. (a) The spin texture of the pristine 3D TI surface is shown 
by the vectors. The circular dotted lines represent the Fermi surfaces 
and the momenta are represented in units of A ^. (b) The corre¬ 
sponding current distribution is shown in the momentum space for 
the injected electron having its spin along the x (red), y (blue) and 
the 2 (cyan) direction, (c) A schematic of the proposed setup to mea¬ 
sure the current asymmetry is shown. The symmetric line (dashed) 
about which the two contacts (shown in blue) are placed makes an 
angle 7 with the x axis as shown. The current A is the total current 
carried by momentum modes which have their azimuthal angles be¬ 
tween 7 and 7 + TT, and current I2 is carried by momentum modes 
whose azimuthal angles lie between 7 -f tt to 7 -f 27 r. The current 
asymmetry is defined as AI — li — E- (d) The current asymmetries 
are plotted as a function of 7 for the three different spin directions of 
the injected electrons. 


IV. SIGNATURES OF SURFACE POTENTIAES IN 
SPIN-POLARIZED TUNNELING 


The current asymmetry measurements carry signatures of 
the spin textures as discussed in Sec. Ill and the spin textures 
can carry signatures of the surface potentials as mentioned in 
Sec. II. Hence, a natural question is can the current asymmetry 
measurements directly characterize the surface potentials and 
this section discusses this question in details. It turns out that 
the form of the spin textures mentioned in Eq.9 leads to an 
easy identification of the functions in Eqs.lO via the current 
asymmetries and it is indeed possible to fully characterize the 
surface potentials from them. 

The three components of the spin texture can be uniquely 
identified by three scans of the current asymmetries over 7, 
each with the spin of the electron from spin polarized source 
along three orthogonal directions. From the form of the spin 
textures in Eq.9, it can be seen that the current asymmetries as 
a function of 7 for the injected electron spin being along x, y. 
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and z are given by 

A&'> = = -(-A’e sin 7 + <¥ 30087 ) (17a) 

= -'S')^/wf(7) = -(-3^c sin 7 + 3^3 cos 7 ) (17b) 

A/(^) = = -(-ZcSin 7 + Z 3 COS 7 ) (17c) 

jT y J jj- 

Hence, in an experiment one needs to fit the functions in Eq.lO 
to the experimentally measured current asymmetries to recon¬ 
struct the surface potentials present. 

In order to understand the effect of the surface potentials on 
the current asymmetries, the classes of surface potentials, one 
which rotate the spin texture out of the plane, and the other 
which rotates the spin texture within the plane are discussed 
separately. Finally how the current asymmetries for a generic 
case with all surface potentials present with different mag¬ 
nitudes can reconstruct the surface potential configuration is 
discussed. 


A. Surface potentials which rotate the spin out of the plane 

The spin textures are rotated out of the plane of the surface 
if the axis of rotation, A lies in the plane (in this case, the x-y 
plane); the surface potentials of the form A^, ^ 0 , Ay ^ 0 , 
and A 2 = 0 fall under this category. The component of the 
in-plane spin texture along A stays the same, however the in¬ 
plane component perpendicular to it gets reduced in magni¬ 
tude. Accordingly, the out of plane component of the texture 
picks up a finite value. 

Owing to the azimuthal symmetry of the Fermi surface, 
without loss of generality we can analyze the case of A^; 7 ^ 
0 , Ay = 0 = A 2 which rotates the spin texture about the 
x-axis and hence does not affect the spin texture along the x- 
direction This is mathematically manifested as A), = 0 and 
<¥3 = 1 for this case. However the functions and 3^3 devi¬ 
ate from their pristine values of -1 and 0 respectively, and so 
do Zc and Zg from 0. This means that the spin textures along 
the y and z direction are affected and this is also evident from 
the example shown in Fig.2a. It is clear that the component of 
the spin along the y direction for every momenta gets reduced 
and the out of plane z component picks up a finite value. 

Since the surface potential does not affect (ct“), AI mea¬ 
sured with the injected spin along the x direction {Os = 
'kI2,(^S = 0 ) has the same behavior with 7 as the pristine 
case. However, as the component of the spin texture along 
the y direction is reduced in magnitude, the value of Aj for 
the injected spin along y direction {6s = 7 r/ 2,())5 = 7r/2) 
at its extrema is also correspondingly reduced as shown in 
Fig.2d. In fact, the extremum value for carries a direct 

information of the magnitude of the surface potential A^, as 
7 r /2 = 23^c/ 7J‘- The out of plane magnetization (cr^) 
also picks up a texture corresponding to the effect on the {a'^) 
due to the surface potential. Hence the current asymmetry 


with injected spin along the z direction also shows an extrema 
at 7 = 7 r /2 whose value is given by 2 Zc/ 7 r. To avoid the sys¬ 
tematic errors in an experiment it is better to look at the ratios 

AI<yy'>\^=^/ 2 /AP\^=o and whose 

experimentally obtained values can be fitted to the curves of 
3^c and Zc shown in Fig. 2c to uniquely determine the magni¬ 
tude of the surface potential A^,. 

At this point, it is important to mention some particular fea¬ 
tures of the behaviors of yc and Zc- Note that in the limit 
of Aj, ^ m, 3 ^c and Zc tend to their values of the pristine 
texture i.e -1 and 0 respectively. This implies that in the limit 
Ax ^ m, its effect on the spin textures is negligible. A^^ = m 
is another special value where 20 a = tt which implies Zc = 0 
(its pristine value), whereas 3^0 = 1 which is negative of its 
pristine value. This leads to the spin texture looking like as if it 
is obtained from a Dresselhaus spin-orbit coupling. A change 
of sign in AJ^^^ as a function of 7 compared to the pristine 
case signals this. For A^, = {y/2 ± l)m, it turns out that 
20 a = TT 12 which means the spin texture lies completely in 
the x-z plane as manifested in 3^c = 0 and Zc = +1. This situ¬ 
ations would reflect itself in the current asymmetries as 
being identically zero for all 7 and |(A/(^^ l 7 =o)| 

taking its maximum value 1 . 

For completeness we also present the results for the case 
of Ay 7 ^ 0 , Ax = 0 = A 2 , which rotates the spin texture 
about the y-axis. As a result the spin texture along the y di¬ 
rection remains unaffected and along the x direction it reduces 
in magnitude. As before the out plane magnetization picks 
up a texture in accordance to the loss of magnetization along 
the x direction. It can be seen from Fig. 3d that the current 
asymmetry profile for the injected spin along the y direction 
is the same as the one for the pristine case which indicates that 
{a^) has the same texture as the pristine case i.e yc = —1 and 
3^3 = 0. However Xg and Zg deviate from their pristine values 
of 1 and 0 respectively, leading to and A/^^) showing 

behaviors different from the pristine case. 

B. Surface potentials which rotate the spin in the plane 

The surface potential of the type Azcr^ ® has an effect 
on the spin texture which is quite different from the previous 
case as it rotates the spin for all the momenta by the angle 
20a about the z axis and hence changes the resultant spin- 
momentum locking angle. Contrary to the previous case, it 
leaves the spin texture completely planar as it does not in¬ 
duce any out of plane magnetization on the surface states of 
the 3D TI. The changed spin-momentum locking angle can be 
directly obtained from the current asymmetries, as for a par¬ 
ticular spin direction of the injected electron, the extrema in 
AI occurs at ^max = (j^s — (^L + TT/2. In the presence of 
the surface potentials the maximas of AI get shifted by angle 
A as shown in Fig.4d. In an experiment the value of A can 
easily be obtained as mentioned and it could be used to di¬ 
rectly obtain the value of the surface potential A^ by fitting A 
to sin -1 ( 3 ^ 3 /y 3 fT 3 f). 

It is also interesting to note that at A^ = (v^ ± l)m, 
20a = TT 12 which means the spin texture becomes completely 
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FIG. 2. (a) The spin texture is shown for = 0.3m, = 0 = IS.^ 

with the vectors showing the in-plane spin textures {a^) and (ct^) 
and the color density plot in the background showing the out of plane 
spin texture (cr^). Note that the magnetization along the y direction is 
reduced in magnitude, (b) For the same A, the current distribution is 
shown as a polar plot for three different spin directions of the injected 
electron. It can be seen that the maxima of the current distribution 
for the case of the injected spin along y direction is reduced. Also, 
due to the spin texture along the « direction, the current distribution 
for the injected spin along the z direction is no longer circular. These 
facts are also reflected in the current asymmetries shown in (d). The 
behavior of the functions and as a function of A^, is shown 
in (c). The inset shows their behavior in the limit of large A^, where 
they tend to their pristine values. 

radial with (cr) = ±r. This will be reflected by the AI as 
function of 7 getting shifted by 7 r/ 2 . The spin-momentum 
locking angle 9l as a function of A 2 is shown in Fig. 4c 


C. Generic case 

In the previous two subsections the effect of the surface po¬ 
tentials which rotate the spin-texture in and out of the plane 
were studied separately. However, in a realistic generic case, 
it is expected that all of them would be present together with 
different magnitudes. In such a case, their effect on the spin 
texture could be a complicated combination of the effects dis¬ 
cussed in Secs. IV A and IV B. It is then imperative to uniquely 
extract the values A^,, Ay, and A 2 from the current asymme¬ 
try measurements that are accessible experimentally. 

It turns out that this is indeed possible via a constrained 
optimization with the constraints being 

AJ(") 1^=0 = ‘2Xs/tt, (18a) 



0 10 7r/4 7t/2 37r/4 tt 

Ay/m j 


FIG. 3. (a) The spin texture is shown for Ay = 0.3m, Aj, = 0 = Az 
with the vectors showing the in-plane spin textures {a^} and (ct^) 
and the color density plot in the background showing the out of plane 
spin texture (cr^). Note that the magnetization along the x direction 
is reduced in magnitude, (b) For the same A, the current distribu¬ 
tion is shown as a polar plot for three different spin directions of 
the injected electron. It can be seen that the maxima of the cur¬ 
rent distribution for the case of the injected spin along x direction 
is correspondingly reduced. Also, due to the spin texture along the 
2 : direction, the current distribution for the injected spin along the 2 
direction is no longer circular. These facts are also reflected in the 
current asymmetries shown in (d). The behavior of the functions Xs 
and Zs as a function of Ay is shown in (c). The inset shows their 
behavior in the limit of large A^ where they tend to their pristine 
values. 

= -2Vc/^; AI^y'> 1^=0 = 2Vs/7r, (18b) 

A/(") |^=,/2 = -2Z,/Tr; A&^ |^=o = (18 c) 

Essentially this optimization is effectively a search for a 
triplet of values for A^,, Ay, and A^ which simultaneously 
satisfy the six constraints in Eq.l 8 , which can be done by sev¬ 
eral numerical techniques. 


V. DISCUSSIONS AND SUMMARY 

One of the crucial requirements for the practical implemen¬ 
tation of our proposal is the ability to manipulate the spin- 
polarization of the tunneling tip. The outstanding progress in 
the field of direct observation of spin textures and magnetic 
order via spin-sensitive scanning probes like spin-polarized 
STM and magnetic exchange force microscopy indeed makes 
our proposal feasible (for a detailed review, see Ref.36). It 
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FIG. 4. (a) The twisted but in-plane spin textures for Az = 

0.3m, Ax. — 0 = Ay is shown, (b) The corresponding current 
distribution is shown as a polar plot. Note that the distribution com¬ 
pared to the pristine case shown in Fig.l is rotated corresponding 
to the twist in the spin textures, (d) The twist in the spin texture 
is reflected in the current asymmetries as the profile is shifted in 7 
by an angle A which gives the new spin-momentum locking angle, 
(c) The spin-momentum locking angle is plotted as a function of the 
surface potential strength A^. Note that for large values of A^ the 
spin-momentum locking angle tends to its pristine value of — 7 r/ 2 . 


was found that the polarization of spin-polarized tips pre¬ 
pared by coating non-magnetic tips (Tungsten (W)) with anti- 
ferromagnetic materials (Chromium (Cr)) can be manipulated 
by varying the thickness of the Cr layer.^®’^^ For instance, W 
tips coated with 25-45 monolayers of Cr are sensitive to the 
out-of-plane polarization where as with ca.lOO monolayers 
of Cr are sensitive to the in-plane polarization.^®'^® This pro¬ 
vides a way of manipulating the polarization of the tip with¬ 
out changing its orientation to its crystalline axis or the crystal 
growth axis of the 3D TI. In our case, with the tip being sen¬ 
sitive to the out-of-plane polarization, the z component of the 
spin texture can be measured. In order to scan the in-plane 
spin texture, the in-plane polarization of the tunneling tip can 
be fixed and one does not need to necessarily rotate it. The 
point to note here is that the tunnel magnetoresistance sig¬ 
nal simply depends on the relative angle between the contact 
orientation and the tip polarization (Eq.l 6 in the manuscript). 
By having multiple pairs of contacts with different orienta¬ 
tions (denoted by the angle 7) on the same sample, the current 
asymmetry can be measured as a function of 7. 

Note that in our case, with the in-plane polarization of tip 
fixed, one could also rotate the orientation of the contacts rel¬ 
ative to the tip polarization, or alternatively, taking advantage 
of the azimuthal Fermi surface, rotate the sample with the con¬ 


tacts fixed so as to effectively change its orientation relative to 
the tip polarization. 

The experimental feasibility of our proposals also relies on 
the ballistic nature of transport on the surface of the 3D TI. 
Hence, the elastic mean free path is a crucial length scale 
in the problem. Since the transport studied in the work is 
spin polarized, the spin relaxation length is a more relevant 
length scale, however in the absence of magnetic impurities it 
is bounded from below by the mean free path®®. The elas¬ 
tic mean free path of 3D TI material Bi 2 Te 3 has been re¬ 
ported to be 235nm^®, so it is expected that a sample with 
a micrometer-by-micrometer surface would be sufficient for 
the experiment and it is indeed a feasible sample size'®. Our 
proposal also relies on the contacts being reflectionless and as 
different amounts of reflections from the two contacts could 
lead to a skewed current asymmetry. However, the presence 
of time reversal symmetry and spin-momentum locking on the 
surfaces of the 3D TI ensure that the contacts are largely re¬ 
flectionless as any reflection would require a spin flip which 
is not possible by non-magnetic contacts. In fact experiments 
with contact resistances as low as ImH have been reported®^. 

In this work, the hexagonally warped regimes of the 
spectrum"" have not been explored, however the derivation 
of the matrix A4 is exactly the same in the presence of the 
warping term. So, the effect of the warping can be very sim¬ 
ply taken into account by taking the appropriate form of \iPq) 
(in Eq. 6 ) to include its effect. In fact, the combination of the 
effects warping and surface potentials can lead to exotic spin 
textures on the surfaces. 

It is also important to mention that though our calculations 
have been done for surfaces of 3D TIs perpendicular to the 
crystal growth axis, they can be straightforwardly be extended 
to surfaces at arbitrary angles (9) to the crystal growth via the 
transformations®®’®® which take a ^ Sg and t ^ Tg where 
Sg and Tg are 0-dependent linear combinations of the spin 
and orbital pseudospin given by 

Tg = {aTy; + /3Ty (g) (Ty , ttTy - g CTy , 

Sg = {aax - g cr^,ay, g a^;}, 

where ^3 = y/ (w|| sin 0 )® + {vz cos 0 )®, a = Vz cos 9/vs, and 
P = U|| sin0/v3. 

It is also worth mentioning that our protocol can be ex¬ 
tended to study weak disorder which varies slowly on the sur¬ 
face of the 3D TI provided the disorder potentials preserve 
time reversal symmetry, however a rigorous treatment of dis¬ 
order on the surface of 3D TIs presents an interesting problem 
in itself. 

Einally, to summarize the main results of the work, a sys¬ 
tematic treatment of the all possible time reversal symmetric 
surface potentials on 3D TIs was presented and their effect on 
the spin-momentum locking of the surface states was studied. 
It was found that these surface potentials can have non-trivial 
effects leading to novel spin textures on the surfaces of 3D TIs. 
Quantum transport experimental protocols were suggested to 
study these spin textures using spin polarized electron injec¬ 
tion into the surfaces of the TI leading to multi-terminal tun¬ 
nel magnetoresistance response. It was shown that these stud¬ 
ies can be used to extract all the information quantitatively 
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about the surface potentials, if present, on the surfaces of 3D 
TI and hence can act as a smoking gun for them. These find¬ 
ings are extremely important from the point of view of the 
spin-momentum locked surface states of 3D TIs being used 
for spintronics applications. Finally, the surface potentials can 
also serve as a platform for engineering novel spin textures on 
the surfaces of 3D TIs, provided a controlled way of manipu¬ 
lating the surface potentials in experiments is engineered. 
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Appendix: Derivation of the matrix A4 


In this appendix the derivation of the matrix A4 is sketched. The derivation involves taking the Hamiltonian with the surface 
potentials included and integrating the Dirac equation across the boundary between the 3D TI and the trivial vacuum to obtain 
the discontinuity in the spinors due to the surface potentials as they are modeled by a Dirac-delta function. As mentioned in the 
main text, the spinor on the trivial vacuum side (z > 0) stays the same as \ijjo) which we denote as |'0o) | 2 =o+ = cg, dg)’^. 

The spinor on the TI side (z < 0) is denoted by \fpA)\z=o- = (^A, &A, ca, c^a)^ 

The Hamiltonian H in Eq.4 in its matrix form is given by 


/ 


0 


ike 


-mg -Vzdz - ils,z5{z) 

a • A c/ \ ike~'^^^v\\-\- 

zdz+l^z5{z) mg ,(A,.-.A„) 5 ( 2 ) 

^ —ike^^^vw— ^ .A f-/ \ 

0 -"^0 -Vzdz+l^z5{z) 

V z(A.+*aJ5(2) 0 

such that when the Dirac equation is integrated across the boundary at z = 0 as 


i(Aj;—iAy)<5(z) 

0 


mg 


r+t 

dz = lim / dz E\'ijj), 

-e 

it gives rise to a system of four linear coupled equations which are given by 

2vzibg - 6a) + i^z{bg + 6a) + « - iAy'j {dg + cIa) = 0, 

2vz(ag - qa) + iAziag + oa) + i - iAy^ (cq + ca) = 0, 

— 2vz{dg — c?a) + iAz{dg + (Ia) — i ^Aa; -1- iAy'j {bg + 6a) = 0, 

- 2vzicg - ca) + iAz{cg + Ca) - i (Aj; -I- iAy'j {ag + ua) = 0. 
These equations can be rearranged to express the entries of |'0a) in terms of the entries of It/jg) as 

= M 



where the matrix A4 is given by 

/_|AI2 I --.2 

1 

M = 


1 

+ |A|2 


V 


TO" — 2iAzm 0 

0 —\A\'^+m?' — 2iAzm 

0 

0 


(A.l) 


(A.2) 

(A.3a) 

(A.3b) 

(A.3c) 

(A.3d) 

(A.4) 


2AyTO — 2iAxm 


2AyTO — 2iAxm 


—2Aym — 2Axim 0 

0 —2AyTO — 2Axim 

— \A\^+m? + 2Azim 0 

0 — I AP-I-TO^-I-2A2ZTOy 


(A.5) 


which can be written in a simpler form given in Eq.7. 
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